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Abstract 

We study fluctuation effects in a two species reaction-diffusion system, 
witli three competing reactions yl + ^—>0, 5 + 5^0, and A + B ^ 
0. Asymptotic density decay rates are calculated for d < 2 using two 
separate methods - the Smoluchowski approximation, and also field the- 
oretic/renormalisation group (RG) techniques. Both approaches predict 
power law decays, with exponents which asymptotically depend only on the 
ratio of diffusion constants, and not on the reaction rates. Furthermore, we 
find that, for d < 2, the Smoluchowski approximation and the RG improved 
tree level give identical exponents. However, whereas the Smoluchowski ap- 
proach cannot easily be improved, we show that the RG provides a systematic 
method for incorporating additional fluctuation effects. We demonstrate this 
advantage by evaluating one loop corrections for the exponents in d < 2, and 
find good agreement with simulations and exact results. 
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1 Introduction 



Over the past decade there has been enormous interest in reaction-diffusion systems 
(see [1-12] and references therein), with particular emphasis on the effects of fluc- 
tuations in low spatial dimensions. Most attention has been paid to reactions of the 
form A + A ^ and A + i? ^ with a variety of different initial/boundary condi- 
tions. At or below an upper critical dimension dc, these systems exhibit fluctuation 
induced anomalous kinetics, and the straightforward application of traditional ap- 
proaches, such as mean field rate equations, breaks down. Attempts to understand 
the role played by fluctuations for d < have involved several techniques, includ- 
ing Smoluchowski type approximations and field theoretic methods [§, [10], |lT 



In this paper we set out to study these fluctuation effects in a system with three 
competing reactions: 

The reactions are irreversible, and we choose homogeneous, though not necessar- 
ily equal, initial densities for the two species at t = 0. Our goal is to calculate 
density decay exponents and amplitudes, taking into account fluctuation effects. 
In pursuit of this aim, we analyse the system using both the Smoluchowski ap- 
proximation and the field theory approach, and we show that the two methods are 
closely related. However, whereas it is unclear how the Smoluchowski approach 
may be improved, the field theory provides a systematic way to obtain successively 
more accurate values for the asymptotic density decay exponents and amplitudes. 
We shall concentrate on situations where one of the two species is greatly in the 
majority (as is almost always the case asymptotically) - so, for example, if the A 
species is predominant, then we can safely neglect the reaction B + B ^ ^. This 
kind of assumption will lead to a considerable simplification in our analysis. 

Previous work on this problem includes use of the Smoluchowski approximation 



0, as well as exact Id results obtained by Derrida et al. |15, 16, for the special 
case of immobile minority particles. Derrida et al. were, in fact, studying a different 
problem, namely the probability that a given spin has never flipped in the zero 
temperature Glauber dynamics of the q-state Potts model in one dimension. By 
solving that model exactly |16|, they showed that this probability decreased 
as a power law: t~^^^ for the Ising (g = 2) case. However, in one dimension, 
the Ising spin flip problem and the decay rate for the immobile impurity in our 
reaction-diffusion system are exactly equivalent problems, and hence this exact 
decay rate also holds in our case. We also mention one other previous result for 



the immobile impurity problem, due to Cardy |jT8[. Using renormalisation group 



methods similar to those employed in this paper, it was shown that the density of 
the minority species decays away as a universal power law: for d < 2, where 
/3 = i + 0(e) and e = 2 - d. 

The case where the majority species is immobile has also been solved (see ||19|| ). 
In this case the decay rate for the minority species is dominated by minority impu- 
rity particles existing in regions where there happen to be very few of the majority 
particles. Since these majority particles are strictly stationary, this situation is 
not describable using a rate equation approach, and it turns out that the minority 
species decays away as exp jf:'^/('^+2)j ^ result which is not accessible by pertur- 
bative methods. 

In this paper, using a field theory formalism and techniques from the renormal- 
isation group, we will obtain decay rates and amplitudes for the general case of 
arbitrary diffusivities - a regime previously only accessible using the Smoluchowski 
approximation. Our basic plan is to map the microscopic dynamics, as described by 
a master equation, onto a quantum field theory. This theory is then renormalised 
(for d < 2), and the couplings (reaction rates) are shown to have 0(e) fixed points, 
whose values depend only on the ratio of the species' diffusion constants. Note 



that this system (with irreversible reactions) is particularly simple in that only the 
couplings (and not the diffusivities) are renormalised. The next step is to group 
together Feynman diagrams which are of the same order in the renormalised cou- 
plings - i.e. diagrams with the same number of loops. These diagrams are then 
evaluated and a Callan-Symanzik equation used to obtain improved asymptotic e 
expansions for the densities. In this fashion, quantities of interest may be system- 
atically calculated by successively including higher order sets of diagrams (with 
more loops) in the perturbative sum. 

One consequence of the theory is that the asymptotic decay rates and ampli- 
tudes for d < d(. will be independent of the reaction rates - a result which is in 
accordance with the Smoluchowski approach. In fact, all physical quantities below 
the upper critical dimension asymptotically depend only on the diffusivities and 
the initial densities, and in this sense they display universality. 

We now present a summary of our results for the density decay rates. In 
what follows we define n^, ns to be the initial density of A, B particles, and 
5 = {Db/Da) < 1 to be the ratio of the diffusion constants. For d < dc = 2, 
ua ^ ub and Uj^^'^Da^ ^ t ^ ti (where ti is a crossover time derived in section 
4), we have (as in |P]): 



For the minority species, we find, from the RG improved tree level approximation 
in the field theory: 



These decay exponents are identical with the Smoluchowski results. Performing 
a strict e expansion on this RG improved tree level result gives an exponent (3 = 
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i + 0(e) for the immobile impurity case (5 = 0). This is in agreement with previous 



RG calculations by Cardy [|I8|. If we now go beyond the tree level calculation by 
including one loop diagrams, then we obtain an improved value for the exponent 
P using an e expansion: 
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is the dilogarithmic function pO[ . This exponent is found to be in good agreement 



with simulations HSIl and exact results iTBI in d = 1. 



However, for 5 < 1, the system crosses over to a second regime where (6) ^ (a). 
This situation is similar to the case where we begin with ub ^ n^. In that regime, 
at times Dst ^ ^b^'^i ^"^^ ^r ns ^ riA, 6^0 and c? < 2, we have: 
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for the majority species. Using the RG improved tree level result for the minority 
species, we obtain: 

(a) ~ E{DBt)-'', (7) 
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The exponent is again in agreement with the Smoluchowski result. If we attempt 
to improve this calculation to one loop accuracy, then we obtain: 
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This exponent is only valid for S quite close to unity, and even in this region it may 
be less accurate than the (non e-expanded) RG improved tree level result given 
above. This point will be discussed further in section 4.2. 

We next give results valid for d = 2, where we find extra logarithmic factors 
multiplying the power law decay rates. Treating first the case (a) 3> (6) and 5 < 1, 
we have, from the RG improved tree level, an initial regime with: 

, , \nt , , 

However, for 5 < 1, the system again crosses over to a second regime where (b) ^ 
(a) . In this second regime the density decay exponents (though not the amplitudes) 
are the same as for the case where we begin with hb ^ tia- In that case we have, 
for 5^0: 



Int 



<«) = O I ( — j I . (13) 

Grossover times for these cases are given in section 4.3. 

We now give a brief description of the layout of this paper. In the next section 
we analyse the system using the mean field/Smoluchowski approach. We then 
set up the necessary formalism for our field theory in section 3, and use it to 
perturbatively calculate values for the density exponents and amplitudes in section 
4. Finally, we give some conclusions and prospects for future work in section 5. 



2 The Mean Field and Smoluchowski Approach 

The simplest description of a reaction-diffusion process is provided by the mean 
field rate equations. For the system we are considering with densities a and b, they 
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take the form: 

da 

It 
db 



— = -2\AAa^ - ^ABab (14) 



-2XBBb^ - XAsab, (15) 

where A^^, ^bb, and A^^ are the reaction rates, and where we impose initial 
conditions of the form a\t=o = ua and b\t=o = ub- In this approach we have 
completely neglected the effects of fluctuations - in other words we have made 
assumptions of the form (ab) oc (a)(6) etc., where the angular brackets denote 
averages over the noise. Below the critical dimension, where fluctuations become 
relevant, this sort of approximation will break down. 

Nevertheless, even at the mean field level, the complete solution set for these 
rate equations is quite complicated. In what follows we shall restrict our analysis 
to the case where 2Xbb < ^ab < 2Xaa- The solution for this particular parameter 
set will be required for our later field theoretic analysis. Following f^, it is easy 
to show (by forming a rate equation for the concentration ratio) that (a/6) as 
t oo. Thus if we begin with initial conditions where ua ub, we can identify 
two distinct regimes - an early time regime where a ^ b and, after a crossover, a 
late time (true asymptotic) regime where b ^ a. Treating the early time regime 
first, we find (after some algebra): 

a ~ i2XAAt)-' (16) 

6 (17) 

{2nAXAAt)^^^/^^^^' ^ ' 

Note that the A particles are decaying away more quickly than the B's, so eventually 
we crossover to a second regime: 



6 ~ {2XBBtY^ 



Ua I, , (Xab - 2Xaa) riA] ^^bb(>^ab-'^>-aa) 



(2nBABBt)^-^s/2Ass \^ {2Xbb - ^ab) ub ^ 



■ (19) 



Alternatively, if we begin with ub ^ n^, then we have a single asymptotic regime: 

b ~ i2\BBt)-^ (20) 

However, if we now wish to extend our results at or below the upper critical 
dimension, we must attempt to include some of the fluctuation effects. The simplest 
way in which this can be done is to employ the Smoluchowski approximation p!3|, P, 



1^. The essential idea of this approach is to relate the effective reaction rates A|fjj 
to the diffusion constants Da, Db- Suppose we want to calculate the reaction rate 
^AB- We begin by choosing a (fixed) A species target "trap", which is surrounded 
by B particles. When a B particle approaches within a distance R of the target, 
a reaction is deemed to have occurred. Consequently, the reaction rate may be 
obtained by solving a diffusion equation with boundary conditions of fixed density 
as r oo, and absorption a.tr = R. The flux of B particles across the d dimensional 
sphere of radius R is then proportional to an effective microscopic reaction rate. If 
we now generalise to the case where both the A and B species are mobile, then we 
find (in dimension d < 2 and in the large time limit): 

A^^i ~ (const.)(Z}^ + DBY^h'/'-'. (22) 

For d = 2 we obtain logarithmic corrections: 

^ HiD, + DB)t) ■ ^2'^ 
The Smoluchowski reaction rates for A^^ and \^bb obtained in a similar fashion. 
Note that above d = 2 the reaction rate approaches a limiting (constant) value, 
and we see that the Smoluchowski approach predicts a critical dimension of dc = 2 
for this system. This is simply related to the reentrancy property of random walks 
in d < 2. It is the inclusion of this effect which accounts for the improvement 
introduced by the Smoluchowski approach. 
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If we now substitute these modified reaction rates into the rate equations, we 
can obtain the Smokichowski improved density exponents. For the case where 
Ub, we find an initial regime with 

a = 0{t-^/^) (24) 
h = 0{t—^^—) . (25) 



Once again, since the A particles are decaying away faster than the B's, we cross 
over to a second regime, where (for < 5 < 1) 

h = 0{t-''l^) (26) 



d ( 1+15-1 \ 



d/2\ 



a = 0\^ '\ ' J j. (27) 

This second set of exponents is the same as for the case where we begin with 
ub ^ riA and 5^0. In this situation no crossover occurs and the exponents 
are valid for all asymptotic times. These exponents can be compared favourably 
with both simulations 0, and exact results [0. For example, the decay rate for 



an immobile minority impurity is given by Smoluchowski to be ~ t"*^'^^"*. This 
compares well with the exact decay rate of t""''^^^. 

Turning to the case d = = 2 and n^i ^ n^, we obtain, for the initial regime: 

a = O i^—j (28) 

6 = 0(^(^y'*'\lnt)mK^)j. (29) 
We again eventually crossover to a second regime, where (for < 6 < 1): 

b = oM (30) 



t 

I 



a = 



(31) 



This second set of exponents is again valid (for all asymptotic times) in the case 
where we begin with ub ua and 5 7^ 0. 
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Note that the Smoluchowski approach can also be employed for d > dc, where 
again we will find (time independent) reaction rates which depend on the diffusion 
constants. However, in the general case, our later field theoretic analysis shows 
that there is no real justification for this procedure. One exception to this occurs 
in the case where we have heterogeneous single species annihilation, as considered 
in ip. In this situation we have only one fundamental reaction process, but dif- 
ferent reaction rates may still arise, for example, by having two or more different 
particle masses (and hence two or more different diffusion constants). In this case it 
is physically reasonable to suppose that the exponents (which are ratios of reaction 
rates) may depend only on the diffusivity ratios, with any other parameters cancel- 
ing out. However, in the general case, where the reaction processes are genuinely 
distinct this will not be the case. 

Overall, we have seen that the Smoluchowski approach is a simple way to in- 
corporate some fluctuation effects into the rate equation approach. Unfortunately, 
it is not at all clear how these methods may be systematically improved. It is for 
this reason that we turn to the main purpose of this paper - the development of an 
alternative fleld theoretic framework. 



3 The Field Theory Approach 

Fluctuation effects in react ion- diffusion systems have previously been successfully 
tackled using techniques borrowed from quantum fleld theory and also from the 
renormalisation group. Examples include studies of the diffusion limited reactions 



A + ^ m and A + B ^ ^ lH, 0]- The flrst step in this analysis is to write 



down a Master Equation, which exactly describes the microscopic time evolution 



of the system. Using methods developed by Doi |^ and Peliti ||2^, this can 



be mapped onto a Schroedinger-like equation, with the introduction of a second 
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quantised Hamiltonian, and then onto a field theory, with an action S. These steps 



have been described in detail elsewhere I^T], ^ ^ |10|, |TT|, and consequently we 



shall simply give the resulting action appropriate for our theory: 

S = J d'^x (^j dt [d{dt - V^)a + b{dt - 6V'^)b + 2XAAaa'^ + \aao^c? (32) 
+2XBBbb'^ + ^Bsb^b'^ + XABO'Cd) + XAsbab + XAsdhab — an a — buB^ ■ 

Here we have defined S = {Db/Da) < 1 and also introduced the response fields d 
and b. In addition time t, together with the reaction rates A{jj} have been rescaled 
to absorb the diffusion constant Da- Averaged quantities are then calculated ac- 
cording to 

= J VaVdVbVbX{t) e"^, (33) 

where 

M = j VaVdVbVbe-^. (34) 
Notice that in the path integral 

j VaVdVbVbe-^ (35) 

integration over the fields a, a and b, b whilst neglecting the quartic terms, leads 
to a recovery of the mean field rate equations. 

Performing power counting on the action S, we can now give the natural canon- 
ical dimensions for the various parameters appearing in the action: 

[a],[6],M,M -A;"^ [a], [&] ~ [A{,,}] ~ A;^-^. (36) 

Notice that the reaction rates become dimensionless in d = 2, which we therefore 
postulate as the upper critical dimension for the system, in agreement with the 
Smoluchowski prediction. 

From the action S, we can see that the propagators for the theory are given by 

r e-'^'^*-*') for t > t' 
Ga-a{k,t-t') = \ (37) 

I for t < t' 
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f e-'='(*-*')'5 for t > t' 
G,^{k,t-t') = \ (38) 

I for t < t'. 

Diagrammatically, we represent Gaa by a thin solid line and Gf^i by a thin dotted 

line. The vertices for the theory are given in figure 1. 
3.1 Renormalisation 

One of the most important features of this theory, as mentioned in the introduction, 
is the relative simplicity of its renormalisation. Examination of the vertices given in 
figure 1 reveals that it is not possible to draw diagrams which dress the propagators. 
Hence the bare propagators are the full propagators for the theory. Consequently, 
the only renormalisation needed involves the reaction rates X{ij}i and in particular 
the diffusion constants (or 5) are not renormalised. 

The temporally extended vertex functions for the reaction rates are given by the 
diagrammatic sums given in figure 2. As is the case in similar theories |TU|, P |, 
these sums may be evaluated exactly, using Laplace transforms: 

= l + A,,CT(e)2)(. + iP)-^/^ ^''^ 

\AB{k, ^) = 1 + AAs2-^/2Cr(e/2)(l + (5)-<^/2(s + ^5/(1 + 5))-^l^ ' ^^^^ 

where C = 2/(87r)'^/^ and s is the Laplace transformed time variable. 

We can now use these vertex functions to define the three dimensionless renor- 
malised and bare couplings, with s = k^, k = as the normalisation point: 

9Ri,,y = '^''>'{ij}{k,s)\s=^2^k=o ^0{,,} = (42) 

Consequently, we can define three /? functions: 

d 

PidRAA) = ^g^dRAA = -^9Raa + ^CT{e/2)gl^^ (43) 
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d 



Pign^s) = '^^dR.B = -^9Ra. + 2-^/^eCr(6/2)(l + 6)-'/'gl^^, 
and three fixed points Pig^^..^) = 0: 

9*R,, = (CT(e/2))-i 
9*Rss = (CT(6/2)r 



(44) 
(45) 



-d/2\-l 



9r, 



CT(6/2; 



1 + 



(46) 
(47) 

(48) 



Finally, we see from (pQ]), (^0]), and (^Tj) tliat tlie expansion of go^^^^^ in powers of 
9R{^j} is given by: 

(49) 



9 R^j^jy 

go^ijy = gRi^jy + — + 

9Rw} 



3.2 Callan-Symanzik Equation 

We now exploit the fact that physical quantities calculated using the field theory 
must be independent of the choice of normalisation point. This leads us to a 
Callan-Symanzik equation: 

•9 ^, , d ^, , d 



"'dg 



Raa 



{a)n = 0. (50) 



However dimensional analysis implies 



d ^ d , d J d , 

K- 2t— + duAT: \- duBT. " 

OK ot on A OUb 



{a)R{t,nA,nB,gR^^^^,5,K) = 0. (51) 



Exactly similar equations hold for Eliminating the terms involving k and 

solving by the method of characteristics, we find: 

{a)fi(t, riA, ub, gRi^^yJ, k) = (K^t)"''/^(a)/j(K"^ fi^(K"2), ?2b(k"^), ^^j^^^j (k"^), 5, k), 

(52) 

with the characteristic equations: 



2t^:r— = —dnA 2t^:r— = —duB 2t — = p{g 



dt 



dt 



dt 



R 



(53) 
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and initial conditions: 

n^(t) = Ha fiBif) = ub (54) 

'dRAAit) = 9Raa gRBsit) = 9Rbb 9RABit) = 9Rab- (55) 
These equations have the exact solutions: 



/ t\ ' f t\ ' 

Mt') = [j;) UBit') = (^-j Ub, (56) 



and 



In the large t limit ^ij^^^j — > 5'/?^ a relationship which will allow us to relate an 
expansion in powers of the renormalised couplings QRi^^^y to an e expansion using 
(|52D. In our later density calculations we will assume that this asymptotic regime 
has been reached. 



3.3 Tree Diagrams 

In order to perform systematic e expansion calculations we now need to identify 
the leading and subleading terms in an expansion in powers of Qo^^^t^ ■ In calculating 
(a) and (6), contributions from tree diagrams are of order 5'0{,j}^{i}''5 integer 
g, and densities = {n^jn^}. However, diagrams with / loops will be of order 
"^^^ addition of loops makes the power g^i^^.^ higher relative to the power 
of the densities - so we conclude that the number of loops gives the order of the 
diagram. 

The lowest order diagrams contributing to (a) and (6) are the tree diagrams 
shown in figure 3. We represent the classical (tree level) density {a)ci by a wavy 
solid line, and {h)ci by a wavy dotted line. These sets of diagrams are equivalent to 
the mean field rate equations, as may be seen by acting on each by their respective 
inverse Green functions. 
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The second tree level quantities appearing in the theory are the response func- 
tions: 



L{k, t2, h) = {a{-k, t2)a{k, ti)) (58) 

M(k, t2, h) = {b(-k, t2)a(k, h)) (59) 

N{kMM) = {h{-kM)KkM)) (60) 

P{k,t2,h) = {a{-k,t2)b{k,h)) (61) 

which we represent diagrammatically by the thick lines shown in figure 4. These 
functions can be evaluated analytically, but only in the limit (a) ^ (b), or (6) ^ 
(a). The details of this calculation are presented in appendix A, where the following 
results are derived (for (a) » (6)): 

N(k, t2, h) = (^^l^f^) exp (-e(t2 - m (63) 

P(A;, t2, ti) = -AABriA Tr~W{ 7"^^ (^2 - *!())) 

(^i + ZAAAnAh) 

M{k,t2,h) = -A^^ns^^-^^^^^^^^^^ 

^ A, (1 + 2A W)^ ■ ^ ^ 

An extra check on vahdity of these response functions is provided by the relations: 

P(0,,0,.^|M M(0.,.0)^^. (67) 

which follow from the definition of the response functions, and from the initial 

condition terms in the action S. It is easy to check that the above response functions 
do indeed satisfy these relations. 
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For the opposite situation where tia (and hence (6) ^ (a)), we could use 

a formahsm similar to the above for the density calculations. However, it is much 
simpler to map this case onto the (a) ^ (6) regime by swapping the labels on the 
A and B particles, and then relabeling: 



We can then obtain the exponents and amplitudes for this second regime with no 
extra work. 

This concludes our discussion of the field theory formalism. The framework we 
have built up allows (in principle) the systematic calculation of fluctuation effects 
in all circumstances. However, it is only in the case where one of the species is 
greatly in the majority where the equations (for the tree level densities and response 
functions) are sufficiently simple for analytic progress to be made. We now turn to 
use of the field theory in calculating the fiuctuation modified densities. 

4 Density Calculations 
4.1 Tree Level 

The first step in using our field theory to include fluctuation effects is to insert the 
mean field (tree level) solution into the Callan-Symanzik solution (^), using the 
results for the running densities/couplings (|56|), (pTf). Since the fixed points for the 
couplings obey '^QRgg < Qr^b < '^Qraa (^^en 5 < 1) it is appropriate to use the 
mean field solutions derived in section 2. For the case where ua ^ ub, this gives: 



ua ^ ub Xaa ^ Abs Da ^ Db- 




(68) 



and 



(6) ~ FiDAt)-^ 



(69) 
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with 

valid for n^^^^D^^ <^t <^ti, where 

ri. D I 

DAh 



(71) 



These modified crossover times are obtained by using the expressions for the run- 
ning couphngs/densities in the mean field crossovers. Notice that the density decay 
exponents derived here are the same as those obtained from the Smoluchowski ap- 
proach. However, as we are performing an e expansion, we are only strictly justified 
in retaining leading order e terms. Consequently we find, for the minority species 
density decay exponent and amplitude: 

Eventually, however, as the A particles are decaying away more quickly than the 
B particles (due to their greater diffusivity when 6 < 1), we crossover to a second 
regime where (6) ^ (a). For < 5 < 1, we have: 

(i-> ~ (^) (DBtr"'^ (73) 
(o) ~ E(DBt)"-, (74) 

with 
where 

, [((l + ^)/2)-^/^-lW \ V ' ^ I J 
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This result is valid for t ^ t2, where 

Ul+5-^)/2y/2 ■ U^J 

Note that for 5 = 1 the first crossover time ti — > oo - in this case the two species 
decay away at the same rate, and so no further crossover occurs. Alternatively if 
5 = 0, then the first regime is left, but the second crossover time ^2 oo. In that 
case the minority species finally decays away in the exponential fashion predicted in 
[0. For the intermediate case where 5 is small, but nonzero, the decay exponent 
for the minority species becomes large in the final regime. The explanation for 
this result lies in the relatively large diffusivity of the minority A species (if is 
large) and/or the increased density amplitude for the majority B particles (if -D_b is 
small). Both these effects will lead to an increased rate of decay for the A species. 

Finally, if the initial conditions are changed such that now ub ^ ^a, with 
5 7^ 0, then we obtain the same results as for the second of the above regimes for 
BbI > u^b''^. with / ^ 1. 

4.2 One Loop Results 

We now describe the one loop improvements to the tree level result. In the regime 
(a) ^ (6), the dominant diagrams will be those where the minimum possible 
number of (6)c« insertions are made. For the majority A species the appropriate 
diagram is shown in figure 5, where there are no (6)ci insertions. This is identical 
to the one loop diagram for yl + A evaluated in f^, which gives, in conjunction 
with the sub leading terms from the tree level: 

In addition, for that subset of diagrams with no (6)ci insertions, the decay expo- 
nent is exact. More details of this calculation, including a demonstration of the 
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cancellation of divergences, can be found in p. 

Turning now to the one loop calculation for the minority species, the appropriate 
diagrams are the three shown in figure 6, each of which contains just one {h)ci 
insertion: 

m -^XAB\\An\nB f d^k rt 



(2AAAnAt)^^s/2^^>i J {2tiY Jo Jo Jt, (1 + 2AAA^At2 
^ — exp [-A;2(t2(l + 5) - 2ti + (1 - (5)t')] 



(1 + 2XAAnAt^ 



{2XAAnAtV^^/^^^^ J {27iY Jo "Jo '{l + 2XAAnAt2y 

xexp[-A;2(i + 5)(t2 -ti)]. (82) 

The detail of the evaluation of these diagrams is rather subtle. Essentially 
we are interested in extracting the most divergent parts of these integrals, which 
will turn out to be pieces of 0{e~^) and 0(e°). However, we must be careful not 
to confuse genuine bare divergences (of 0(e^^) which must be removed by the 
renormalisation of the theory), with logarithmic pieces, which we must retain. The 
divergences arise in diagrams (i) and (iii) as the difference in time t2 — ti between 
the beginning and end of the loops tends to zero (in d = 2). After the process of 
renormalisation we find corrections of the form: 

1 + (constant)eln((constant)t°'/^) + O(e^). (83) 

If this series is identified as the expansion of an exponential, then we find that our 
one loop diagrams (together with subleading components from the tree level) have 
provided 0(e) corrections to the exponents. 

Diagrams (i) and (iii) are relatively straightforward to evaluate. The k and ti 
integrals are elementary, and the final t2 integrals can be done by parts to extract 
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the necessary most divergent pieces (up to 0(e°)). The second diagram of figure 6 
is more comphcated, and we perform its evaluation in appendix B - although we 
are only able to extract the logarithmic piece of 0{t^'^^'^^'^^'^^ t'^^'^lnt). There will 
be corrections to this of 0{t~^^^^'^^'^^ t"^/^) (contributing to a modified amplitude) 
which we have been unable to calculate. We find asymptotically: 

8.(2A^^n^t)A.W^A.. [ + t (H^XAAUAt) - 1) ln(8.) 

15t^/2 3 ^ rt 



-f /2 H2XAAnAt) - ^ h'-"'^^ ln(l + 2XAAnAt2)dt2 + 0(e) j (84) 



4 



(85) 



+ 0{e)]f/^H2XAAnAt) 



fi+s , (i +V)^ , , 

- / -^In 1 + f 

J-i 

(m) 2Ar(2ALnU--/--- I 1 + ' H2XAAnAt) ln(4.(l + ^)) 

(86) 

-f/\H2XAAnAt) - 1) - ln(l + 2XAAnAt2)dt2 + 0(e)) . 

To one loop accuracy we can make the replacement: A{ij} = k'^qq^..^ — * K'^g^^^.y 
These results must now be combined with the subleading terms from the tree level. 
Using (^, we find 

(M = !^ ( 1 3R,s 



(2AA^n^i)^AW2Ax^ (2«:^(7K^^nAt)^«-^^/'^«^^ V S^^I^^a 

2 \ 

^^^'^ ln(2«:^^?«,,n^t) + |^ln(2«:^^?^,,n^t) + 0(^?^) . (87) 



If we now insert explicit e expanded values for the fixed points then we dis- 

cover that the bare divergences cancel between (@), (|86D, and (|^). With insertion 
into the Callan-Symanzik solution (|52D , we also find that the pieces we have left as 
integrals in (i) and (iii) (which are 0(t'^/^(lnt)^)) also mutually cancel. Eventually 
we find: 

(6) ~ (const.)r ^ 1 + ^ 



4 8 \ Vl + 5 
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{l + vf 
dv 



1 



ln(l + v) 



ln((const.)t'^/') + 0{ 



where we have neglected 0(e) pieces which, aside from the prefactor, are time 
independent. These terms contribute only to the density amplitude. We now 
evaluate the integral in (p8|), using 



T+l ln(l + v] 



dv 



du 



u 



1 In-u [T+a In-u 

-du + / -du 

u — 1 Ji u — 1 



6 11 + 5 



^9) 



where /{x} is the dilogarithm function |^0[. The other parts of the integral are 
elementary. The next step is to e expand the RG improved tree level result: 



( 1 + 5 ^ 



d/2 



x2/ \ 2 ^ ^ 
Then, exponentiating the e expansion in (|88|), we find (6) = 0(t~^), where 



(90) 



/5 



1 



2 \ 2 4 



2 In 



^(5^-1) (l + (l + 5) 



/ 



1 + 5 



(91) 



f3 is plotted as a function of 6 for e = 1 (ci = 1) in figure 9. For the case where 
5 = 1, we recover the decay rate (6) = 0{t~'^^'^). This is to be expected, as when 
(5 = 1 we are effectively again dealing with a single species reaction-diffusion system 
(at least for d < 2). In that case the density decay exponent is known to all orders 
in perturbation theory and is in agreement with our result. For the case where 
5 = and d = 1, the decay exponent is also known exactly to be (6) = 0(t^^'^'^^) 
[0. This can be compared with our result, where we find 



16 4 



64 



0.39 



(5 = 0). 



(92) 



Consequently, this answer is a modest improvement over the Smoluchowski result 
derived in section 2, and also in |^. 
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For the case ^ (and hence (6) ^ (a)), we could follow the same route 
as described above, by evaluating the one loop diagrams shown in figures 7 and 8. 
However, as we mentioned in the last section we can much more easily obtain these 
corrections by swapping the labels on the A and B particles, and then relabeling: 



A 



D 



B- 



Following this procedure, the majority species amplitude/exponent can be found 
by taking Da ^ Db m equation (|79|) : 

^ 1 21n87r-5 



+ 



0(6) (Dst) 



-d/2 



(93) 



y47re IGvr 

We can obtain the one loop minority species exponent by substituting 6 in 
equation (pi]): 

(a) = 0(t-"), (94) 



where 



a 



'1 + 6-^^ 



1 



^ 2 



3 , 1 + d- 

- + In 

2 I 2 



l + 21n 



1+6- 



1) i + d + r^: 



/ 



(95) 



1 + 5-1 

Notice, however, that in forming the one loop corrections for the minority species 
exponent, we have had to expand the RG improved tree level result: 



1-- i + i„ 



0(r) 



(96) 



The error arising from this expansion will become large as 5 becomes small. Eventu- 
ally this inaccuracy will cause the exponent to reach a maximum and then decrease 
as 5 is further reduced - behaviour which is clearly unphysical. In order to reduce 
the error, and to ensure that the expansion in equation ( PU] ) is qualitatively cor- 
rect, we need to retain the O(e^) terms. Hence the one loop exponent in equation 
(p5D should be treated with some caution - terms of order O(e^) will probably be 



22 



required for precise results. Consequently the (non e expanded) RG improved tree 
level result given in the last section may be more accurate in this regime. In figure 
10 we have plotted the one loop exponent a as a function of 6, for d = 1 (e = 1), 
in the region 0.7 < S < 1, where the exponent is still increasing for decreasing 6. 

In principle, calculations can also be made for the case with ua ^ ub, but 
where we have crossed over to the regime (b) ^ (a) (for < 5 < 1 and times 
t ^ t2). However, a rigorous evaluation of the one loop diagrams is now much 
more difficult, as the functional forms for the densities and response functions 
will change over time. Nevertheless, since the above corrections to the exponents 
come from asymptotic logarithmic terms, it is plausible to suppose that the new 
exponent corrections will be dominated by contributions from the final asymptotic 
regime. If this is indeed the case, then the one loop exponents (though not the 
amplitudes) will be unchanged from the previous results (equations (^) to (P^D). 
This calculation will, however, suffer from the same problem as described above. 

4.3 d = dc 

For the case d = dc = 2 we expect logarithmic corrections to the decay exponents, 
as the reaction rates are marginal parameters at the critical dimension. We 
can find the running couplings from the characteristic equation (^) by taking the 
limit e — > in equations (||), (H), and (H): 



^--(-")- i+,Jcr^Mn(.^t) -(^^"^-^)" 

where we have taken the asymptotic limits. Corrections to the asymptotic run- 
ning couplings will be an order (lnt)~^ smaller, and consequently these asymptotic 
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expressions will only be correct at very large times. Hence our expressions for 
the densities will only be valid when both this condition, and the crossover time 
constraints given below, are satisfied. In what follows we shall assume the validity 
of the first of these two conditions. Notice that the asymptotic running couplings 
are still ordered Igng^ < gn^g < '^Qraa for 5 < 1, so we can use the mean field 
solutions derived in section 2 as the basis for the RG improved tree level exponents 
and amplitudes. Making use of the Callan-Symanzik solution (^) and the above 
running couplings, we find for (a) ^ {h): 

, , Int , , 

where G = exp ( ( 1 ~ (^+'^)gflAA j j jg a non-universal amplitude correction. 



Note that the next order terms for the minority species are suppressed by a factor 
of only (In In t)/ (In t). Using our expressions for the running couplings/densities 
in the mean field crossovers, we find that these expressions are valid for times 
D^^n^^ Int < t < Ti, where 

2 



(D^Ti/lnTi) ^ [^^-^^ j ■ (102) 

For the case 6 < 1 the system will eventually enter a second regime, where now the 
B species will be in the majority. We have (for 5 7^ 0): 

(a) (104) 

with 



9Rbb \ 9Rab ) 

This is valid for times when t ^ T2, where 



H = exp[ 1 - ^^^^^ K= 1 + — ] . (105) 

\9Rbb \ 9Rab J J \ Ub, 



6- 
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Alternatively, if we begin with ub riA, then for 5 7^ and (Dsi/lni) ^ n^^, 
we have the same results as for the second of the above cases, with K ^ 1. Inter- 
estingly, the logarithmic corrections we have derived in this section using the RG 
approach differ slightly from the Smoluchowski results given in section 2. 

5 Conclusion 

In this paper we have made a comparison of two methods for treating fluctuation 
effects in a reaction- diffusion system. We have found that the Smoluchowski and 
field theory approaches are rather similar - the Smoluchowski approximation, for 
d < 2, giving the same exponents as the renormalisation group improved tree level 
in the field theory. In addition, we have gone on to calculate the field theoretic 
one loop corrections, which have yielded improved values for the exponents. The 
advantage of the field theory is that it provides a systematic way to calculate these 
corrections - a procedure which is lacking in the Smoluchowski approach. Further- 
more the use of renormalisation group techniques has demonstrated universality in 
the asymptotic amplitudes and exponents, in that, for d < 2, they only depend on 
the diffusivities and the initial densities, and not on the reaction rates. 

The theory we have developed in this paper can easily be extended to slightly 
different situations. Consider first an annihilation/coagulation reaction-diffusion 
system, where the following reactions occur: 

A + A^A B + B^B A + B^$. 

The Smoluchowski approach differs from before only in the absence of factors of 2 
in the rate equation terms describing the same species reactions. Consequently, if 
we begin with ^ then the minority species will decay as 

b = O {d < 2). (107) 
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On the other hand, the field theory description lacks only the factors of 2 in the 
action (^). If this difference is followed through then the decay exponent in the 
RG improved tree level is seen to be the same as in the Smoluchowski approach. 
However, this difference of a factor of 2 has a major effect on the response functions 
(where this factor appears as a power), and as a result the new one loop corrections 
will be different from those calculated in section 4.2. These results should be 
compared with the exact solution |2j, ^ for the minority species decay rate 



b = 0(t '^), where: 



^ = 2cos-nr/(l + ^)) - 



Note that in this case, although the Smoluchowski answer is qualitatively correct, 
it deviates considerably from the exact answer. Hence we can see that application 
of the Smoluchowski approach does not always lead to accurate exponents. 

Another possible extension is to consider reaction-diffusion systems with more 
than two species of particle. For example, examining a three species system, we 
could have the reactions: 

Analysis of this situation is very similar to before, and we merely remark that in 
the appropriate asymptotic regimes the Smoluchowski and RG improved tree level 
exponents (consisting of ratios of diffusion constants) are once again identical. 
Hence the convergence between the Smoluchowski exponents and those obtained 
from the RG improved tree level is fairly robust, and is not simply confined to the 
two species systems we have previously been considering. A further possibility is 
to analyse the case where we have a continuous distribution of diffusivities, but 
with only a single reaction channel. This has been studied from the Smoluchowski 
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point of view by Krapivsky et al. 0, and it would be interesting to extend our RG 
methods to include this situation. 

Our theory could also be employed to consider clustered immobile reactants - a 
generalisation of the 5 = case included in our calculations. This situation has been 



analysed by Ben-Naim [|12[, using the Smoluchowski approach, where the dimension 
of the cluster di was found to substantially affect the kinetics. Specifically, for 
CO dimensionality d — (i/ < 2 (in a space of dimension d) a finite fraction of the 
impurities was found to survive, whereas ioi d — di >2 the clusters decayed away 
indefinitely. The formalism we have presented in this paper could be adapted to 
study this clustered impurity problem, where calculations could be made without 
reliance on the Smoluchowski approach. 
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6 Appendix A: Response Functions 

Obtaining an exact analytic expression for the response functions is, in general, 
very hard. Suppose we define the "trunk" to be the line of propagators onto which 
the density lines are attached, as shown at the bottom of figure 4. Difficulties arise 
from diagrams where the "trunk" changes from one propagator into the other, and 
then back again, as shown in last of the diagrams for the L response function in 
figure 4. If diagrams of this type are initially excluded then progress can be made. 
Consider first the two subseries shown in figure 11, for the functions ^(/c, ^i) and 
9{k, t2, ti), where diagrams of the above kind have been excluded. These series can 
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be summed exactly (using the same technique as described in 0), giving: 

t2, h) = exp [-k\t2 - h)) exp (^- j\^^AAa + \ABh)dt^ (109) 

0{k, t2, h) = exp [-k\t2 - ti)6) exp (^- J^AXBsb + XABa)dt^ . (110) 

The full response functions are now given by the diagrammatic equations shown 
in figure 12, where all possible diagrams are included. Written out explicitly these 
give: 

L{k, t2, h) = ^{k, t2, h) - Xab T ^{k, ts, T)a{T)M{k, r, h)dT (111) 

M{k, t2, h) = -Xab T e{k, t2, r)b{T)L{k, r, ti)dT (112) 
Jtl 

N{k, t2, h) = e{k, t2, h) -Xab 0{k, t2, T)b{T)P{k, r, h)dT (113) 



P{k, t2, h) = -Xab I ' i{k, ta, r)a(r)iV(A;, r, t^)dT. (114) 
Jtl 

In general this set of coupled integral equations is intractable - however we can 
make progress in the limit where (a) ^ (6), or (6) ^ (a). Considering the case 
where (a) ^ (6), the dominant contributions to the response functions come from 
diagrams with the minimum possible number of {h)ci density line insertions. Ac- 
cordingly, we can now truncate the full diagrammatic equations, as shown in figure 
13. Notice that to this order L, N , and P contain no {h)ci density insertions, 
whereas M must contain one such insertion. In this approximation we can now 
perform the integrals inside the ^ and 9 functions, using the appropriate mean field 
density: 

^„ ^^ABTIA „ , /l + 2AAAnAt2\'""^''"" 



/ {AXBBh + XABOjdt ^ / , of^^ dt = In 

Jt, Jt, 1 + ZXAAUAt 



Itl Jtl 1 + 2XAAnAt V 1 + '^XAAnxti ) 

(116) 
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and therefore 

2 

1 + 2XAAnAt2 



^k, h, h) = ( ] + l^^^^'^l' 1 exp - h)) (117) 



Using these expressions, it is now straightforward to derive the response functions 
given in equations (p^]), (§1), and 



7 Appendix B: A One Loop Integral 

For the case where (a) ^ (6) the hardest of the three diagrams of figure 6 to 
evaluate is (ii) - see equation (|8TD. We shall evaluate it first in = 2, and then 
deduce its form in = 2 — e. Notice that the extra integration resulting from the 
{b)ci insertion in the loop ensures that this diagram is not divergent. Taking the 
asymptotic part of the ti and t' pieces, we find: 



2\AA{2\AAnAt)^^^/^^^^ J (27r)2 7o 'Jo 'Jt, {l + 2\AAnAt2yt'^ 

X exp (-A;^(t2(l + 5)- 2ti + (1 - 5)t')). (119) 
The k and t' integrals are elementary, giving 



-^AB^B /■* dt2{2XAAnAY f'^ ^2 (^20) 

M\AA{2\AAnAtY^^I^^^^ Jo (1 + 2XAAnAt2)^ 7o ' ' ^ ' 



X 



(t2(l + 5)-2ti) 



1 1 



1-S , ( 2U 
+ . .^ T-TTT In 



(t2(l + 5) -2ti)2 \{i + 5)h)) 
Although the first part of the ti integral is straightforward, the second piece in- 
volving the logarithm is more difficult. However, if we make the transformation 
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we find: 



l-S 



[ dt.rr^^-^r^. In = i(l - 5^)t, dv^-^ HI + V), 

Jo '(^2(1 + 5) -2ti)2 \{1 + S)t2j 8^ ''J-i ^ 

(122) 

where all time dependency has been removed from the integral limits. The final ^2 
integral is then easy to perform, and we end up with: 

'l-5\ 



In 



1 + 6J 

(123) 



\n(2XAAnAt). 



-pd.ii±^ln(l + .) 

J-i 

However, we now need to extend this analysis to determine the behaviour of the 
integral in d = 2 — e. If we take the asymptotic part of all the pieces inside the 
integral, and perform power counting, we find that it should scale as t-'^^As/sAAA 
However, this procedure is not strictly valid, as in moving to the asymptotic version 
a false ^2 = divergence is created. Nevertheless, the integral is dominated by 
contributions from late times where arguments based on power counting should be 
valid. Hence in d = 2 — e we find: 



-XabI^B (r, l/r2 



327rA^A(2AAAnAt)^^«/'^^^ V 2 



1 + Sj 

(124) 

(1 + ^)^ 



+ 0{e)]f/Hn{2XAAnAt). 



-pd.ii±^ln(l + .) 
J-i 

Further subleading corrections (in time), which we have not calculated, will lack the 
logarithm factor, and so will contribute to the amplitude for the minority species 
density. 
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• • 

n^5(t) nB5(t) 

Figure 1: Vertices for the field theory. 




Figure 2: The temporally extended vertex functions (a) XAA{k,s), (b) XBB{k,s), 
and (c) XAsik, s). 
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<a>. 





<b>. 



Figure 3: Tree level diagrams for the densities (a) and (b). 



L = 



P = 



DDDD 



M= DDDD 



N= DDDDDDDD 



where, for example 



Figure 4: The response functions. 
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Figure 5: One loop diagram for (a) (when (a) » (b)). 





Figure 6: One loop diagrams for (b) (when (a) ^ (6)). 



□ □ □ □ □ □ 



Figure 7: One loop diagram for (b) (when (b) ^ (a)). 
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(iii) 



: □ □ □ 



□ 



Figure 8: One loop diagrams for (a) (when (6) ^ (a)). 
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Figure 9: The one loop density decay exponent f3 for the minority B species {{b) = 
0{t~^)) as a function of 5 . 
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b) nil 




Figure 11: The diagrammatic equations for (a) ^ and (b) 9. 
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Figure 12: The full diagrammatic equations satisfied by the response functions. 
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DDDD III 

Figure 13: The truncated diagrammatic equations for the response functions, vahd 
for (a) > (b), or (b) > (a). 
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